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Abstract 



We obtain self-consistent macroscopic equations describing interlayer Joseph- 
. son effect and intralayer disequilibrium in one-dimensional array of Josephson 

o . 

coupled layers. We show that "nonequilibrium coupling" can lead to effective 



spatial and time synchronization and formation of coherent dynamic resis- 
tive state (collective Josephson effect) in Nb-A10x.-Nb stacked junctions and 
HTSC (intrinsic Josephson effect). We propose it to be the origin of collective 



^ , switching phenomena observed in PbBiSrCaCuO. 



It is well known that dc and ac Josephson effects take place in layered superconductors 
when external current flows in direction perpendicular to the layers. Two types of layered 
structures are currently experimentally investigated: artificial Nb-A10a;-Nb stacked junctions 
[|l|] and natural layered high-Tc superconductors 0. This systems show similar dynamic 
behavior and can be considered on common background. Electromagnetic theory of 
layered structures with thick enough superconducting layers is developed in Ref. [Q, but in 
a case of thin layers some other effects are to be considered, especially disequilibrium which 
can be essential if layer thickness do is smaller than characteristic length of disequilibrium 
relaxation Ie This criteria is obviously fulfilled for "quantum" structures with layers 

of atomic thickness (HTSC). For "classical" artificial structures 



Here is the inelastic electron-phonon scattering time, D = Ivp/^ is the electron diffusion 
coefficient, = (87r^T^^r)/(7C(3)) is the energy gap, / is the free pass, vp is the Fermi 
velocity, r = (Tc — T)/Tc. Thus < Ie can be fulfilled at least for r ^ 1. Last experiments 
^ manifest nonequilibrium effects in layered structures. 

The origin of nonequilibrium Josephson effect is well established |]5|-^. First of all Joseph- 
son junction is a source of disequihbrium in nonstationary state due to injection of quasi- 
particles. This results in generation of nonzero covariant potential $ = + {h/2e){d6/dt) 
in superconducting layers, where is the electrostatic potential and 9 is the phase of super- 
conducting condensate ($ = in equilibrium state). In nonequilibrium regime an ordinary 
Josephson relation (dip/dt) = {2e/h)V between Josephson phase difference (p = 02 — 9i and 
voltage V = — (f)2 is violated. Instead we have 



Thus disequilibrium modifies interlayer Josephson effect and self-consistent description is 
necessary [0,^. 

The simplest way to obtain self-consistent equations is Ginzburg-Landau approach. In 
this paper we consider thin layer limit do <^ Ie, Xl and temperatures close to Tc, thus the 
discrete order parameter "^nix, y) for n-th layer can be used in a case of weak disequilibrium 
~ const). Static equations of Ginzburg- Landau type for are obtained by Lowrence 
and Doniach (see also Ref. [p!0| -[r3|). Here we use dynamic generalization of their model in 
one-dimensional case \E'„(x,y,t) = \I'n(^) (three dimensional variant can be derived straight- 
forward). Thus we consider the structure with spatial dimensions less than characteristic 
magnetic length and without external magnetic field, so that Josephson junctions are effec- 
tively zero- dimensional. One-dimensional state is also possible in larger systems but stability 
analysis is necessary. Interlayer interaction is assumed to be small (only two neighbor layers 
interact). 




(2) 



Model time-dependent Ginzburg-Landau equation of Lowrence-Doniach form may be 
obtained from relaxation principle 

7(/.| + 2.e0.)v^. = -^^. (3) 

Here we introduce electrostatic potential of n-tli layer Free energy is (we can 

choose the gauge A=0 in one- dimensional case) 



(4) 



where a and b are GL constants, d is in general some dimensional parameter depending 
on a and b definition and coinciding as a rule with spatial period of a structure {d ^ do 
for "classical" S-I-S structures and interlayer distance for "quantum"), r is dimensionless 
parameter of interlayer interaction. The microscopic expression for r is rather distinct in 
"quantum" and "classical" cases. But from phenomenological point of view it is more 
appropriate to express r from equation used below 



Tic^d 



(5) 



where Aoy and ^o|| are parallel to the layers London and coherence lengths at zero temperature 
and Jo = {Aea^dr) / (hb) is Josephson critical current density at zero temperature. In both 
cases if r <^ 1 then only two neighbor layers interact. 



From we obtain equation for dimensionless order parameter ip^ = "^n/ycer/b 

~ ^ ^T*^") + ~ IV'nH^n - e(2V'n " " ^n-l) = 0. (6) 

e = r/r is coupling parameter, it characterizes a ratio of Josephson energy to energy of 



superconducting condensate. In the weak coupling limit e <^ 1, \ipn\ ~ 1 introducing 

h d9n 
2e~dt 



'4'n = e^^", $„ = 0n + 7^-7—5 = d-n+i — we obtain from (||) the expression for covariant 



potential 



2e 

tcLj^^n = e(sin - sin Lpn-i) , (7) 



where tci = to/r, to = {j%)l{ad) is characteristic relaxation time and can be found to be 
(at l^/'l = const) 115| 







where Tq = j2/7CiS)nn''Tct,. Note, that this expression follows from time-dependent 



Ginzburg-Landau equation derived in Ref. ||T^ from microscopic theory in the so-called 



local equilibrium approximation, and this equation coincides with our model equation in 
weak disequilibrium and "discrete" case. 

Expression for superconducting current can be obtained from variation principle (in 3D 
form) and full current between the n-th and n+l-th layers is sum of supercurrent, normal 
current taken in simplest ohmic form and capacity current due to polarization of interlayer 
space 



Jn,n+l = Jc(V'n+lC " C+lV'n) " ('^"+^^ '^"^ _ Q ^{(f)n+l - 0n) = COHSt, (9) 



here Jc = Jqt, Rn,n+i, C are the dimensional parameters. We neglect self charge of supercon- 
ducting layers (which is in general small due to quasineutrality) and thus one-dimensional 
current is conserved (possible peculiarities in dynamics of "quantum" structures with do 



compared with Debye length will be considered in other paper, see also Ref. IT^). We intro- 
duce random resistivity -R„,„+i = Ro/Sn,n+i to account the effects of nonidentical junction 
parameters in simplest way. Using (|^) which is direct consequence of definition we obtain 
from (P) 



nd'^iPn ^r diPn . 1 2e , ^ , , ^2e d . ^ ^ . , , 

^p" ^ + '^"^ + - u;;'-{<l>^+, - $„) - a;-2--($„+i - $„) = j, (10) 



where Up = {2eJc)/ (hC), Uc = {2eRoJc) /h, j = J/Jc- As a result we have a set of coupled 
equations (0), (p!0| ) describing Josephson effect and disequilibrium in self-consistent manner. 
This equations are in fact more general than underlying time-dependent Lowrence-Doniach 
equation (^). Eq. ([^) demonstrates well known result $ oc dw}s and can be derived at 
any temperature. Eq. (0) is a sequence of basic nonequilibrium Josephson relation (H) 



and expression for current between layers (^). Thus we think that proposed in this paper 
macroscopic model is a good basis for describing nonequilibrium Josephson effect in layered 
superconductors. 

Equations (|^, (0) are equivalent to a set of coupled Josephson equations for with 
additional interaction terms 

jW\ = ^?(2smv9„ - smcpn+i - smy^^.i) + fJr] {2 cos ipn-^ - cosipn+i—^, cosv9„_i 

(3 = uo'l/uj'i, t' = ujct, 7] is the novel parameter of nonequilibrium coupling 



V = = = 1 • (13) 

Let us estimate r] for some real structures. First of all Fg is usually large as (Td/TY, Td is 
Debye temperature 

h fTD\' ^ fTD\^ I 10=^ Low T, SC VI ~ lO^! 

~ — I I To ~ ( TfT 

" 10 HTSC Fg-lOO 



T\T J " \T, 



' c 



And for rj we obtain 

ST^r Fo , . 

V = -T ^- (14) 

We find that t]^/^ ~ 1 for HTSC at r ~ 10~^ c^eo ~ lO^^^ad/sec and for Nb - AlO^ - Nb 
structures at do ~ lOOA, R^S ~ 10^^ Q cm^ {R^ - resistivity, S - area of Josephson junc- 
tion). Thus we see that weak disequilibrium can lead to strong coupling between interlayer 
Josephson junctions and to collective dynamic behavior in a case of weak coupling between 
layers (e ^ 1). In a case of nonidentical junctions this lead to synchronization. Simulations 
show that at 1] > 1 and jSrj < 1 coherent dynamic state with ipn{t) = ^p{t) is single and 
stable. Full investigation of this problem is a subject of future work. 

Presented theory is able to explain recent experiments ||^ on resistivity (at current 
perpendicular to the layers) of BiSrCaCuO and PbBiSrCaCuO superconductors. In Bi- 
compound many hysteresises on voltage - current characteristic are observed due to one by 



one switching of interlayer junctions and in PbBi-compound all junctions are switched to re- 
sistive state simultaneously. The estimation of parameter rj for this samples gives rj ~ 0.1 -j-l 
for the first of them and ~ 0.5 -^ 5 for the second one. Thus collective switching may be 
induced by nonequilibrium coupling. 
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